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复杂 非 均 质地 层 的 等 值 线 生成 算法 研究 
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摘 要 : 等 值 线 图 是 复杂 油气 藏 地 质 研究 中 的 一 类 重要 表示 形式 。 针 对 地 层 数据 特点 ， 研 究 和 分 析 了 复杂 非 均 质地 层 
数据 等 值 线 绘制 方法 的 不 足 ， 提 出 一 种 新 的 分 区 方法 ， 实 现 了 带 断 层 的 数据 集合 划分 。 基 于 Delaunay 三 角 剖 分 算法 ， 
并 以 几何 多 边 形 为 锥 形 ， 建 立 了 带 约 束 的 三 角 网 格 模型 。 采 用 约束 Delaunay 三 角 网 的 拓扑 结构 ， 选 用 改进 的 反 距离 加 
权 插 值 算 法 ， OA 并 结合 实际 应 用 要 求 在 VC++ 与 OpenSceneGraph 的 开发 环境 下 ， 
现 了 非 均 质 带 断层 的 地 层 等 值 线 的 绘制 。 通 过 对 国内 某 采 油 厂 的 实际 油气 储 集 地 层 数 据 等 值 线 图 的 绘制 和 可 视 化 ， 证 
明了 方法 的 可 行 性 和 准确 性 
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Abstract: Contour map is an important type of representation in the study of complex reservoirs. Aimed at the characteristics 


of formation data, this paper firstly analyzed the Shortcomings of the contour drawing method and proposed a new partition 


method, which realizes the partition of data sets with faults. Then it established a triangular mesh model based on the Delaunay 


triangulation algorithm and geometric polygons. Using the topology of constrained Delaunay triangulation and improved inverse 


二 _ distance weighted interpolation algorithm, this paper solved the problem of elevation values on both sides of the faults. In the 


development environment of VC++ and OpenSceneGraph, this paper realized the contour drawing of the stratum with 
heterogeneity fault. Through the drawing and visualization of the contour map of the actual oil and gas reservoir data, it proves 
the feasibility and accuracy of the method. 
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(triangulated irregular network) 生成 算法 中 ，Delaunay 三 角 剂 
分 算法 在 地 形 拟 合 方面 较 普遍 ， 常 被 用 于 TIN 的 生成 。 当 断层 
等 值 线 图 是 地 质 研究 工作 中 的 一 种 重要 的 表示 方法 HY, 为 。 参与 TIN 的 生成 时 ， 应 将 断层 视 为 约束 条 件 ， 则 应 考虑 带 限定 

IX. 


确定 石油 富 集 区 的 位 置 、 形 态 、 石 油 的 储量 以 及 钻井 作业 等 担 ” ”约束 条 件 的 Delaunay 三 角 网 ， 即 约束 Delaunay 三 角 网 
供 了 重要 依据 ， 因 此 地 层 等 值 线 图 的 绘制 与 可 视 化 具有 非常 重 (constrained delaunay triangulation, CDT ) 59。 

要 的 价值 。 由 于 地 壳 运 动 过 程 中 产生 的 作用 力 ， 使 地 层 发 生 断 前 ， 生 成 约束 Delaunay 三 角 网 的 算法 有 很 多 ， 大 致 上 分 
裂 ， 沿 着 断裂 面 两 侧 发 生 相 对 位 移 ， 即 断层 。 断 层 使 其 两 侧 数 为 两 步 : 首先 ， 生 成 无 约束 的 Delaunay 三 角 网 ， 其 次 ， 将 断层 


据 不 再 连续 , 给 地 层 等 值 线 的 绘制 及 地 层 构造 可 视 化 带 来 困难 。 ”约束 线段 嵌入 到 无 约束 的 三 角 网 中 ， 并 对 三 角 网 进行 调整 。 有 

人 主 分 布 在 断层 附近 ， 准 确 的 等 值 线 图 。 关 约 束 线段 嵌入 的 算法 比较 多 [ 导 ， 文 献 [7] 中 ， 将 除 断 层 之 外 
可 以 清晰 地 反映 复杂 地 质 构 造 ， 使 专业 人 员 能 够 对 地 层 构 造作 ”的 数据 先进 行 剖 分 ， 随 后 再 把 断层 嵌入 ， 这 样 ， 断 层 线段 跟 三 
出 正确 的 分 析 ， 从 而 带 来 重大 的 经 济 效益 NN。 在 TIN 和 角 网 的 位 置 情况 会 有 多 种 ， 对 每 种 情况 分 别 进行 不 同 的 嵌入 方 
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共 边 的 两 个 三 角形 ， 


来 比较 耗 时 ， 刘 少 兵 四 提 到 的 复 连 通 区 域 三 
法 首先 要 找到 与 约束 线段 相交 或 靠近 约束 线段 的 三 角形 ， 并 把 
个 单 连通 区 域 ,然后 求 取 该 区 域 的 外 边界 ， 


这 些 三 角形 合并 成 


法 ， 算 法 较 复 杂 ， 不 易 操 作 ;， 黄 唱 晶 外 提出 的 “插入 -交换 ” 算 
法 适合 于 任意 区 域 的 约束 齐 分 ， 但 在 断层 影响 域内 三 角形 ， 每 
就 要 判断 一 次 是 否 为 凹凸 多 边 形 ， 


执行 起 
剖 分 算法 ， 该 方 


将 单 连 通 区 域 与 外 边界 组 成 复 连 通 区 域 ， 这 样 就 转换 成 了 简单 
多 边 形 的 三 角 章 分 ， 该 方法 解决 了 约束 线段 嵌入 问题 ， 但 效率 


张 : 


中 5， 付 金华 局 提 出 的 等 值 线 生 成 算法 ， 解 决 了 单条 


不 高 
EA 


断 


线 的 葡 入 问题 ,但 并 没有 对 断层 多 边 


进行 处 理 ， 孙 劲 光 


[9 提出 当 约 束 条 件 为 折线 约束 或 者 


线 约束 时 ， 可 以 将 约束 转 


换 为 线性 单元 集合 , 然后 对 这 些 线性 单元 集合 进行 Delaunay 剂 
分 ， 但 在 转换 过 程 中 ， 要 保证 转换 后 的 线性 单元 集合 与 原来 约 
束 条 件 拟 合 , 在 剖 分 过 程 中 , 若 线段 不 在 Delaunay 三 角 网 格 中 ， 


许多 文献 0419 是 通过 局 部 变换 边 或 者 边界 
分 ， 以 恢复 约束 线段 在 网 格 中 的 存在 。 
上 的 点 而 不 是 原来 约束 条 件 中 的 点 ， 这 种 
元 集合 对 原来 约束 条 件 的 逼近 。 

厂 地 层 离散 数据 为 基础 ， 结 合 


取 点 进行 加 点 旨 
法 加 入 的 是 线性 单元 
方法 不 利于 逼近 线性 让 


Im 


分 算法 在 直线 上 
这 种 算 


针对 上 述 问题 ， 本 文 以 采 


弓 诺 


值 线 生 成 算法 。 


字 拱 各 


盖 


维 可 视 化 的 应 用 特点 ， 提 出 一 种 针对 非 
该 方法 首先 将 原始 地 层 数据 进行 分 区 ， 
区 后 的 数据 分 别 进行 无 约束 的 三 角 剖 分 ， 然后 通过 文献 7 


匀 质 带 断 层 的 


区 


中 提 到 的 约束 线段 嵌入 方法 ， 将 断层 多 边 形 嵌 入 到 三 角 网 中 


通过 判断 断层 多 边 形 的 内 儿 
当 


边界 类 型 ， 对 构成 的 三 角 网 进行 适 
当 的 裁剪 与 优化 ， 最 后 进行 等 值 线 追 踪 。 该 方法 不 仅 实现 了 带 


断层 的 等 值 线 的 绘制 ， 而 
1 ， 带 断层 约束 的 层面 建 模 
地 层 数据 分 区 


满足 


1.1 


对 于 给 定 的 油 藏 地 层 不 规则 离散 点 集合 v 以 及 断层 数据 集 


了 实际 应 用 需求 。 


合 ，y={X,7Zhi=1L23.2， 其 中 X,7 为 离散 点 的 坐标 ，z 为 


离散 点 的 高 程 值 ， = {x,y),j=1 


坐标 。 其 分 区 思想 如 下 : 


2 -m3 其 中 X,Y 为 断层 点 的 


a) 以 离散 点 数据 集合 为 例 , 设 定 分 区 号 wwm ,num 从 1 开始 ， 


任 取 离 散 点 集合 中 的 一 点 为 p1， 
离散 点 数据 集合 中 取出 


王 意 一 点 


将 其 分 区 号 设 为 1， 从 剩 下 的 
为 P2， 假定 一 条 以 点 p1 与 点 


p2 为 端点 的 线段 ， 如 果 该 线段 与 所 
区 , 将 点 P2 分 区 号 设 为 p1 的 
， 并 将 P2 从 原始 离散 数据 集合 中 删除 。 如 果 点 P2 与 该 


那么 点 D2 跟 点 p1 同 属于 一 个 分 
分 


的 断层 数据 线 都 不 相交 ， 


区 
分 区 号 下 的 pn 存在 相交 的 情况 ， 


和 分 情况 进行 判断 : 


(a) 将 点 P2 先 不 做 处 理 ， 依 次 遍历 离散 点 数据 集合 中 剩 下 


的 点 ， 若 剩 下 的 点 中 至 少 存在 一 点 D3, 且 以 点 p3、 


点 D1 为 端 


点 的 线段 与 所 有 断层 数据 线 都 不 相交 ， 那 么 将 p3 的 分 区 号 设 


为 1 的 分 区 号 ， 修 改 分 区 号 完成 后 ， 判 断 以 点 p2、 点 p3 为 端 


断层 线 是 否 相 交 ， 


点 的 线段 与 所 


和 孙 歧 峰 ， 等 : 复杂 非 均 质地 层 的 等 值 线 


删除 ， 否 则 ，P2 跟 pl 不 属于 同一 个 分 区 ; 
(b) 若 遍历 完 所 有 的 离散 数据 集合 , 不 存在 点 p3， 使 得 以 点 


点 P1 为 端点 的 线段 与 所 


断层 数据 线 都 不 相交 ， 那 么 分 


p3、 


区 号 num+l] ; 


b) 依次 循环 吉 历 ， 直 全 


= 


取 遍 所 有 的 离散 数据 点 。 


c) 以 断层 数据 为 例 。 任 取 断 层 数据 集合 中 的 一 点 1 ， 从 


num =1 开 


台 ， 若 该 分 区 内 至 少 存在 一 点 p4 ， 


使 得 以 点 p44 、 


点 万 为 端点 的 线段 与 所 有 的 断 


入 到 这 一 分 区 内 ， 并 将 其 


导线 都 不 相交 ， 那 么 将 点 1 加 


从 原始 数组 中 删除 。 


d) 依次 循环 过 历 ， 


直到 取 遍 所 有 的 断层 数据 点 。 


如 图 1、2 所 示 ， 黑 色 


圆 点 代表 离散 数据 点 ， 红 色 代 表 断 层 


数据 点 ， 图 


1 为 分 区 示意 图 ， 图 2 为 分 区 完成 图 。 


. . ee ee ee ee 。 
@ © ee @ 。 © 。 
记忆 站 ® P4 
图 1 离散 点 示意 图 


一 第 一 分 区 4 一 。 一 > 第 二 分 区 < 一 


图 2 分 区 完成 图 


1.2 ”构造 带 断 层 约束 的 CDT 网 格 


1.2.1 构建 无 约束 的 DT 网 格 
本 文采 用 逐 点 插入 法 按照 


上 述 分 区 ， 构 建 每 个 分 区 内 的 离 


散 数据 点 与 断层 点 组 成 的 无 约束 DT 


参考 文献 0%!91。 
1.2.2 断层 约束 线段 入 网 


网 格 ， 其 具体 思想 步骤 可 


选取 断层 上 相 邻 两 点 组 成 的 线段 ， 假 设 要 幅 入 的 线段 为 


f1f2，f1、j2 在 原始 三 角 网 中 


立 置 如 图 3 Ca) 所 示 , 与 f1f2 


相交 的 三 角形 组 成 的 区 域 为 11f2 的 影响 区 域 。 给 定 两 个 数组 
来 存放 三 角形 顶点 断层 入 网 的 详细 实现 过 程 如 下 : 


right 、~ left 


a) 确定 影响 区 域 。 从 三 


有 


链表 中 查找 与 约束 线段 1f2 


相交 的 三 角形 ， 保 留 三 角形 的 顶点 ， 在 三 角形 链表 中 删除 该 三 


角形 ; 
b) 判 


Ey 


如 果 不 相交 ， 就 将 P2 的 分 


区 号 设 为 与 p1 相 同 的 分 区 号 ， 


将 p2 从 原始 离散 数据 集合 中 


点 与 


直线 的 位 


及 上 


新 a) 中 存储 的 顶点 与 J1f2 线段 的 位 
关系 判定 法 。 若 点 在 ff2 的 右 


关系 。 这 里 
则 ， 把 该 


点 存储 到 vignt 数组 中 , 若 在 f1f2 的 左 侧 ， 把 点 存放 在 1gf 数组 


中 ， 若 点 在 f1f2 Ey 把 该 点 同时 存放 到 right ~、 left 数组 中 ; 
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c) 循环 遍历 三 角 
相交 的 三 角形 ; 

d) 删除 数组 yignt 、 ef 中 的 重复 元 素 ， 按 顺序 排列 ， 并 在 
数组 末尾 添加 头 元 素 ， 组 成 pz 、 ppR 两 个 多 边 形 ， 如 图 3 (b) 
所 示 ， 红 线 代 表 约 束 线段 。 


链表 ,执行 a) b)， 直 到 


i 


不 存在 与 f1f2 


fl 全 PR fl 


PL 


(a) 约 束 线段 影响 区 域 (b) 影 响 多 边 形 (0) 影响 多 边 形 的 三 角 训 分 
图 3 约束 线段 入 网 过 程 


1.3 ”影响 多 边 形 的 三 角 剖 分 

因 内 外 边界 约束 的 三 角 剖 分 在 有 限 元 分 析 、 曲 面 重 建 等 领 
域 具有 广泛 的 应 用 且 效 率 较 高 ， 故 而 本 文 使 用 内 外 边界 约束 的 
算法 ， 对 多 边 形 进行 前 分 E9q。 当 断层 约束 嵌入 完成 后 ,， 对 pr 、 
PR 分 别 进行 三 角 剖 分 ， 若 pL 、 PR 为 凸 多 边 形 ， 不 做 处 理 ， 若 


为 凸 多 边 形 ， 要 将 pz 、 pr 外 部 多 余 三 角形 删除 ， 并 进行 优化 
调整 。 因 此 处 讨论 断层 多 边 形 的 约束 算法 ， 当 完成 上 述 断 层 线 
段 杏 入 ， 影 响 多 边 形 的 三 角 剖 分 后 ， 若 断层 多 边 形 内 还 存在 多 
余 的 三 角形 ， 需 将 断层 多 边 形 视 为 边界 ， 判 断 断 层 多 边 形 的 边 
界 类 型 ， 若 为 内 边界 ， 将 断层 多 边 形 内 部 的 三 角形 删除 ， 若 为 
外 边界 ， 需 将 断层 多 边 形 外 部 多 余 的 三 角形 删除 。 删 除 多 余 三 
角形 的 算法 步骤 如 下 : 
a) 判定 断层 多 边 形 的 边界 类 型 。 


b) 计算 三 角形 的 内 心 。 因 三 角形 内 心 始 终 存在 于 三 角形 内 
部 ， 故 而 ， 判 断 三 角形 的 内 心 与 断层 多 边 形 的 位 置 关 系 ， 若 断 
层 多 边 形 为 内 边界 ， 那 么 删除 三 角形 内 心 在 断层 多 边 形 内 的 三 
角形 ; 若 断 层 多 边 形 为 外 边界 ， 那 么 删除 三 角形 内 心 在 断层 多 
边 形 外 的 三 角形 ; 

(a) 三 角形 内 心 的 计算 公 
标 为 A(xl,yl) 、 


Ww 


\ 式 为 : 假定 三 角形 的 三 个 项 点 的 坐 


C(x3,y3) ， 


B(x2, y2) 、 WB 
内 心 为 及 (xX,y) 。 则 内 心 点 用 (x,y) 的 坐标 公式 由 式 (1)~(3) 


计算 得 到 ; 


AB=c， 


a= N(x2—x3) +(y2—y3) (1) 


式 《3)， 类 似 可 以 推出 p、c 的 值 。 


x=(a*xl+b*x2+c*x3)/(atb+ce) (2) 
y=(a*yl+b*y2+c*y3)/(atb+c) (3) 
(Db) 点 与 多 边 形 位 置 关 系 判 定 方法 。 判断 点 与 多 边 形 的 位 置 


关系 ， 常 用 的 方法 为 射线 法 。 即 从 要 判定 的 点 引出 一 条 射线 ， 
统计 该 射线 与 多 边 形 的 交点 个 数 。 若 有 奇数 个 交点 ， 说 明 点 在 
多 边 形 中 ， 和 否则 ， 点 在 多 边 形 外 。 需 注意 的 是 ， 当 检测 的 点 是 
多 边 形 的 顶点 或 者 是 多 边 形 边 上 的 点 时 ， 在 此 处 算 作 点 在 多 边 
形 内 。 如 图 4 为 点 与 多 边 形 的 位 置 关 系 示意 图 (红色 代表 要 判 


乍 期 刊 


仙人 
孙 歧 峰 ， 年， 旬 示 和 5 FE 由 


定 的 点 ， pl 代表 判定 点 与 多 边 


点 重合 )。 


点 pl 在 多 边 形 上 
| 点 p2 在 多 边 形 内 


p3 


点 p3 在 多 边 形 外 


对 


4 点 与 多 边 形 位 置 关 系 示意 图 


c) 循环 遍历 a) b)， 直 到 所 有 的 断层 多 边 形 处 理 完 毕 。 

以 上 述 数 据 为 例 ,图 5 为 断层 线段 岁入 完成 后 ,形成 的 CDT 
网 格 , 图 6 为 删除 多 余 三 角形 后 形成 的 CDT 网 格 (红色 代表 断 
民 数 据点 ， 黑 色 为 离散 数据 点 )。 


图 5 嵌入 断层 后 的 CDT 网 格 


图 6 删除 多 余 三 


形 后 的 CDT 网 格 


2 ” 带 断 层 的 等 值 线 生成 算法 


2.1 反 距 离 加 权 插 值 算法 
前 ， 经 野外 勘探 获得 的 断层 数据 信息 都 是 不 完整 的 ， 大 
部 分 断层 缺失 或 者 没有 高 程 值 信息 ， 所 以 在 绘制 带 断 层 的 等 值 
线 图 时 ， 需 计算 断层 两 侧 的 高 程 值 。 用 于 等 值 线 绘制 时 网 格 插 
值 的 方法 有 很 多 ， 如 : 反 距离 加 权 法 、 克 里 金 插 值 法 、 最 邻近 
点 法 等 。 本 文采 用 反 距离 加 权 插 值 算法 ， 该 方法 比较 简单 、 容 
易 实现 ， 且 有 效 地 解决 了 计算 断层 高 程 值 的 问题 。 
反 距 离 加 权 插 值 (Inverse Distance Weighting) 又 叫 最 小 二 
乘 距 离 加 权 插 值 ， 它 的 插值 点 的 值 只 是 周围 邻近 点 距离 倒数 的 
函数 PU。 对 于 插值 点 x 的 估计 值 为 z ， 其 一 般 形式 的 公式 为 
zo) =D wz /Dw (4) 


其 中 : 
w (0) =1/ d(x,x)" (5) 


d= Vx + -yy (6) 

在 式 (4) ~ (6) 中 ，x 为 临近 点 的 横 坐 标 ，z 为 临近 点 的 

高 程 值 ，w 为 用 于 插值 的 样本 点 的 个 数 ,( 此 处 N=6 )，4a 为 临 
近 点 x 到 插值 点 x 的 距离 ,1 为 正 整数 寡 次 (一 般 0.5 到 3 的 
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值 可 获得 最 合理 的 结果 ， 此 处 m=2 )，w 为 权重 ， 距 离 插值 点 
越 近 ， 邻 近 点 的 权重 就 越 大 。 邻 近 点 的 选取 ， 通 过 计算 离散 数 
据点 到 插值 点 的 距离 ， 通 过 设 定 距离 长 度 ， 查 找 距 离 插 值 点 最 
近 的 6 个 点 ， 从 而 得 到 邻近 的 离散 点 。 


2.2 等 值 线 追踪 


进行 等 值 线 追踪 时 ， 需 注意 等 值 线 是 否 到 达 断 层 线 ， 若 已 


成 后 的 CDT 网 格 进行 结构 化 处 理 


经 追踪 到 断层 线 ， 应 立即 停止 ，1 
现 穿越 断层 的 情况 ， 在 进行 等 值 线 追 踪 2 


j 不 是 穿越 断层 线 。 为 避免 出 
前 ， 需 对 约束 前 分 完 
， 将 断层 线 所 在 的 三 角形 视 


为 边界 三 角形 ， 当 等 值 线 追 踪 到 边界 三 角形 时 ， 停 止 追 踪 。 追 


踪 算 法 如 下 : 


孙 歧 峰 ， 等 : 


-1520 
-1544 一 一 一 983 
Sw 52 一 435 


双 


中 10 带 断 层 的 等 值 线 


在 断层 线条 数 为 6， 等 值 线 为 11 的 情况 下 ， 对 不 同 数据 级 


别 的 数据 进行 测试 ， 时 间 对 比 结果 如 表 1、 表 2 所 示 。 


a) 遍 历 所 有 的 三 角形 集合 , 查找 边界 三 角形 及 断层 三 角形 ， 
并 将 该 三 角形 中 的 边界 边 及 断层 边 存储 在 边界 边 链表 中 ， 


b) 以 某 等 值 线 数据 为 例 , 从 边界 边 链表 中 选取 一 条 等 值 边 
作为 起 始 边 。 查 找 以 该 边 为 边 的 三 角形 及 其 相 邻 三 角形 ， 并 记 


录 穿 过 的 所 有 三 角形 ， 直 到 


查找 到 的 等 值 边 的 最 后 一 条 是 边界 


Ke 


边 ， 非 闭合 的 等 值 线 追 踪 完 


c) 闭合 的 追踪 与 非 闭合 类 似 , 但 需 注 意 的 是 ， 其 追踪 结束 


条 件 为 查找 到 的 等 值 边 的 最 后 一 条 为 起 始 边 。 


3 ”算法 实现 实例 


本 文 以 油 藏 地 层 离散 数据 为 基础 ， 实 验 环境 为 Visual 
Studio2013+OpenSceneGraph， 实 现 了 非 均 质 带 断 层 的 等 值 线 绘 


格 ; 图 9 为 无 约束 上 


验证 明 该 方法 效果 


-16r1593 


制 功 能 。 图 7 为 无 约束 DT 网 格 , 图 8 为 带 断 层 约束 的 CDT 网 
的 等 值 线 图 ， 图 10 为 带 断 层 的 等 值 线 图 。 实 
良好 。 


一 一 一 一 二 一 


图 7 无 约束 DT 网 格 


一 EN 
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带 断 层 约束 的 CDT 网 格 
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图 9 无 约束 等 值 线 图 


表 1 数据 集合 约束 剖 分 耗 时 
实验 组 别 1 组 2 组 3 组 4 组 
离散 点 数 200 1000 5000 10000 
本 文 算 法 耗 时 /ms 18 120 482 921 


表 2 传统 等 值 线 生成 与 本 文 等 值 线 生 成 算法 耗 时 对 比 


实验 组 别 1 组 2 组 3 组 4 组 
离散 点 数 200 1000 5000 10000 
传统 算法 耗 时 /ms 25 116 580 1160 
本 文 算法 耗 时 /ms 15 89 445 890 
4 ”结束 语 
本 文 在 前 人 研究 的 基础 之 上 ， 对 带 断 层 的 复杂 地 层 等 值 线 
绘制 问题 提出 了 一 个 新 的 方法 。 该 方法 首先 将 地 层 离散 数据 点 


及 断层 数据 点 划分 在 不 同 区 块 中 ， 简 化 了 数据 集 ， 通 过 断层 约 
束 线段 的 嵌入 、 影 响 多 边 形 的 三 角 剖 分 以 及 断层 数据 插值 等 操 
作 ， 完 整 的 实现 了 非 均 质 带 断层 的 等 值 线 的 可 视 化 。 该 方法 为 


有 效 地 进行 地 层 建 模 分 析 与 应 用 葛 定 了 坚实 的 基础 。 
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